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EXECUTIVE SUMMARY 

Future space exploration and commercialization will require more efficient heat re- 
jection systems. In order to accommodate the required heat transfer rates, such systems 
must employ advanced heat transfer techniques. Forced two-phase flow boiling heat trans- 
fer with enhancements falls in this category. However, moderate to high quality two-phase 
systems tend to require higher (compared to single-phase systems) pressure losses. Pre- 
vious investigators have studied a variety of two- phase flow (with F!reon-ll) regimes in 
uniformly-heated (or cooled) flow channels. It has been demonstrated that various degrees 
of heat transfer enhancement are possible depending on the: (1) enhancement device used, 
(2) predominating flow regime, (3) geometry (e.g., orientation, heat flux distribution), and 
(4) specific localized thermophysical processes. 

Since many space applications will involve single-side (top- side) heating, it is essential 
to extend the literature to assess any enhancement of such heating alone and with other 
enhancement techniques. Experiments, using freon-11, were used primarily to examine top- 
side heating with and without enhancements. The objective of this work was to expand 
the existing freon- 11 flow loop capabilities to include both saturated and subcooled flow 
boiling experiments. The experiments include: (1) measurements of local (axial) heat 
transfer coefficients, (2) evaluation of the effectiveness of combined (e.g., combined helical 
fins and twisted tapes) and single heat transfer enhancements, (3) assessment of the effects 
on heat transfer of single-side heating, and (4) development of heat transfer data reduction 
techniques. 

Although most flow boiling experiments documented in the technical literature deal 
with uniform heating of the coolant channel walls, actual cold-plate components will always 
be heated from one side. In addition, the coolant channels of these components are likely 
to have a horizontal orientation. However, the magnitude of the gravitational vector will 
be usually small or nonexistant. Although stratification will not occur in actual cold plate 
channels due to gravity, some form of stratification may occur due to: (1) circumferential 
surface tension gradients, and (2) localized or single-side heat flux. 

Thermal systems designers, for the most part, are using data from uniformly heated 
experiments to design components which are subjected to single-side heat fluxes. Under 
certain flow conditions such heat flux distributions can significantly affect the heat transfer 
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and hence the ability of the component to transport thermal energy. The present work has 
increased the heat transfer data for single-side heated channels with enhancement devices. 

This final report is divided into two major parts: (1) Part I: Multidimensional Wall 
Temperature Measurement and Heat Transfer Enhancement for Top-Heated Horizontal 
Channels with Flow Boiling, and (2) Part II: Improved Analytical Heat Transfer Data 
Reduction for a Single-Side Heated Coolant Channel. 

Part I summarizes over forty experiments which involve both single-phase convection 
and flow boiling in a horizontal channel heated externally from the top side. The following 
experimental parametric effects were included: (1) inside wall enhancement [smooth, spiral 
fins with small pitch, spiral fins with large pitch (LP) combined enhancements which 
include LP and a twisted tape], (2) channel inside diameter, (3) inlet subcooling, and 
circumferential heat flux distribution [uniform and top-side heating]. Part I also has an 
appendix, where the experimental uncertainty analyis was summarized. The resulting 
uncertainty in the mean heat transfer coefficient was ±14.6 W/m?K. 

Part II contains parametric dimensionless curves with parameters such as the coolant 
channel radius ratio, the Biot number, and the circumferential coordinate. These curves 
directly relate the measured circumferentially averaged temperature to an equivalent mean 
heat transfer coefficient. For the case of non-uniform circumferential heat transfer distri- 
bution on the outside of the coolant channel, most designers desire to know the equivalent 
non-uniform heat flux distribution on the inside of the coolant channel. This inside sur- 
face heat flux distribution has been obtained for the case of constant external heat flux on 
the top of the coolant channel with the bottom half insulated. The inside wall heat flux 
distribution is presented in terms of dimensionless radial temperature gradients. 

This report only summarizes many of the activities of the students and the principal 
investigator. More details can be found in students individual reports and theses. This 
work has supported four (4) graduate students and five (5) undergraduate students. How- 
ever, many more students have been exposed indirectly to this project through technical 
presentations to student groups and tours of the laboratory. As a result, more students 
have considered and/or made plans to pursue a graduate education or improve their aca- 
demic performance. 

Finally, the continual support of this work by NASA and funded research from 
other agencies have assisted in the creation of the Thermal Science Research Cen- 
ter (TSRC) at Prairie View A&M University. This Center is interdisciplinary in that 
professors from Engineering, Physics, Chemistry, Mathematics, Computer Science and En- 
gineering Technology will focus their vast backgrounds on thermal science related research 
and design. This includes heat transfer, thermodynamics, thermal control, fluid mechan- 
ics, mass transport, thermo-chemical processes, micro-heat transfer, two-phase flows, and 
non-intrusive optical diagnostics. 
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PART I: 


MULTI-DIMENSIONAL WALL 
TEMPERATURE MEASUREMENTS 
AND HEAT TRANSFER ENHANCEMENT 
FOR TOP-HEATED HORIZONTAL CHANNELS 

WITH FLOW BOILING 
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Abstract 

Two-dimensional (circumferential and axial) wall temperature distributions have been 
measured for top-heated coolant channels with different internal geometries which include 
smooth walls, spiral fins, and both twisted tape and spiral fins. Freon- 11 was the working 
fluid. The flow regimes studied include: (1) single phase, (2) subcooled flow boiling, and 
(3) stratified flow boiling. The inside diameter of all test sections was near 1.0 cm. Cir- 
cumferentially averaged heat transfer coefficients at several axial locations were obtained 
for selected coolant channels for a mass velocity of 210 kg/m 2 s , 0.19 MPa (absolute) exit 
pressure, and 20.8 ° C inlet subcooling. Overall (averaged over the entire channel) heat 
transfer coefficients were compared for the above channel geometries. This comparison 
showed that the channel with large pitch spiral fins had higher heat transfer coefficients 
at all power levels. However, the results appear to indicate that if the twist ratio (ratio of 
the twisted tape period to the inside diameter) was decreased, the configuration employing 

1 Honeywell Endowed Professor of Engineering and Director of the Thermal Science 
Research Center 

2 Harrison Division, General Motors Corp.; 200 Upper Mountain Rd; A&cE, Bldg 6; 
Lockport, NY 14094 

3 Westinghouse Savannah River Company; P. O. Box 616; Bldg 705-lc; Aiken, SC 29802 
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both fins and a twisted tape would have had greater enhancements. 


NOMENCLATURE 


Surface area, m 2 
Thermal hydraulic diameter, m 
Eckert number, t ~ — r 

Ptiq c p\* * $at) 

Circumferentially and axially averaged heat transfer coefficient (HTC), W/m 2 K 
Axially-distributed but circumferentially-averaged HTC, W/m 2 K 
HTC due to natural convection (see Figure 4), W/m 2 K 
Specific latent heat of vaporization, kJ/kg 
Enthalpy subcooling of the fluid, kJ/kg 
Jakob number, — * ut> ** 

Thermal conductivity, W/mK 
Peclet number, Re Pr 
P v Net power generation, W 

q c Heat flux due to natural convection from outside of test section (see Figure 4), W/m 2 

q R Heat flux due to radiation from the outside of the test section, W/m 2 

r Radial coordinate for the data reduction model (see Figure 4), m 

Tf Bulk temperature of the flowing fluid (see Figure 2), °C 

T m (4>,Z ) Local measured outside wall temperature of the test section (see Figures 3 and 4), °C 
T w {<}>, Z ) Outside wall temperature of the test section = T m {<j>, Z), °C 
T S at Saturation temperature (316 K at 0.19 MPa for Freon-11), °C 


A s 

D 

Ec 

h 

hjn 

hoo 

ifg 

Ja* 

k 

Pe 
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Too Ambient temperature, °C 

We* Modified Weber number, ^ 

Z Axial coordinate for the heated portion of the test section (see Figure 3), cm 

Subscripts 

A , B, and C Denotes Domains A, B, or C (see Figure 4) 
i =1 through 7 for each axial location (see Figure 3) 

Greek 

<f> 

7T 

Pg 
Pliq 


Circumferential coordinate; see Figures 2, 3, and 6 (”Phi”) 

Half of a full rotation or 180 degrees; in some figures, 7r is also referred to as "Pi.” 
Density of gaseous phase of fluid, kg/m 3 
Density of liquid phase of fluid, kg/m 3 
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INTRODUCTION 


Space commercialization will require efficient heat transfer systems. The future suc- 
cess of many efforts will be based on our understanding of the behavior of two-phase flow 
boiling in both the space (zero-g or reduced-g) and earth environments. This was empha- 
sized in the Workshop on Two-Phase Fluid Behavior in a Space Environment, sponsored 
by NASA (Swanson, et al., 1989). Flow boiling heat transfer offers an enhancement al- 
ternative to forced-primed and capillarity heat management systems. Essential to better 
understanding factors affecting flow boiling in heated tubes is the following effects: (1) 
non-uniform heat flux distribution, (2) local (axial and circumferential) distributions of 
the heat transfer coefficient, (3) resulting pressure drop and pumping power requirements, 
(4) single and double enhancement devices, (5) the relative advantages of saturated and 
subcooled flow boiling regimes, (6) flow channel aspect ratio effects, (7) the relative effects 
of heat transfer enhancement techniques, and (8) correlations for mean and local heat 
transfer, and pressure drop. In addition to being applicable to several gravitational levels, 
future research efforts must also include basic phenomena such as: (1) orientation (e.g., 
vertical flow and bottom-heated flow channels) and Marangoni effects, (2) other working 
fluids such as ammonia, (3) flow stability, (4) binary fluids, and (5) identification of the 
threshold inertia (Froude number) beyond which gravity effects would be negligible. For 
example, threshold inertia determination is necessary to identify when orientation and/or 
Marangoni effects become important. Although it is not apparent, the development of 
improved data reduction models is also essential to the accurate representation and inter- 
pretation of the heat transfer data. 
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This work will assist the development of fundamentally-based heat transfer corre- 
lations which include effects of: (1) enhancement device configurations for fluids other 
than air (Camovos, 1977), (2) basic flow parameters which are fluid independent, and (3) 
complex heat flux distributions. 

This paper examines, experimentally, two-dimensional wall temperature variations for 
highly turbulent horizontal channel flows, which axe heated from the top. Freon-11 was 
the working fluid. The temperature measurement were used to obtain axial variations in 
the circumferentially averaged heat transfer coefficient {h m ). The over-all heat transfer 
coefficient (h) was determined and compared for four different internal channel enhance- 
ment configurations (see Table 1): (1) smooth wall, (2) spiral fins with a small pitch (SP, 
6.52 fins per cm), (3) spiral fins with a larger pitch (LP, 4.0 fins per cm), and (4) doubly 
enhanced spiral fins with both large pitch fins and a twisted tape. 


EXPERIMENTAL INVESTIGATION 

In this section, brief descriptions Eire given of the freon-11 flow loop, test section, and 
the data reduction procedure. 
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Flow Loop 


The flow loop (see Figure 1) was a closed system, which operates between 0.1 MPa 
( T, a t — 24 °C) and 1.3 MPa ( T sa t = 124®C). The total power generation was 2.6 kW, and 
the maximum mass velocity was 1.0 Mg/rrPs. Under special circumstanses, the loop can 
be operated at a pressure near 0.04 MPa ( T sa t = 0°C). The loop has two reservoirs (0.25 
m 3 each). The reservoirs have separate heat exchanger jackets for secondary temperature 
control. The flow loop tubing was 2.5 cm (inside diameter) stainless steel. The flow 
loop consists of: (1) the freon loop, (2) air lines for pneumatically controlling valves and 
pulsation damper, and (3) vacuum lines for system evacuation. 

The freon loop was designed to study both staurated and subcooled flow boiling 
regimes. After the loop was evaculated, it was filled with freon until the pressure was 
slightly above the atmospheric pressure. After the charging process, the fluid was circulated 
through the loop at the desired operating conditions. Bled valves were used to purge the 
loop and transducers of any gases, and the flow conditions were reestablished. During a 
given test, the heater tape power output was adjusted to a given level. The flow rate of the 
isothermal bath (50%-50% mixture of ethylene glycol and water) was adjusted until the 
test section inlet temperature was at the desired level. A steady-state was then allowed 
to occur. After all the desired flow conditions were again verified, the test section’s axial 
and circumferential wall temperature measurements were recorded along with all the flow 
conditions. 

Referring to Figure 1, the freon flowed form reservoir 1 through a filter to the main 
pump ( positive displacement), which required a net positive suction pressure of at least 


9 



0.02 MPa. After leaving the pump, the fluid passed near a pulsation damper, which reduced 
the pressure an d flow oscillations. When the damper was used, a pneumatically controlled 
metering valve was used to stabilized the flow. After exiting the metering valve, the fluid 
passed through a turbine flowmeter and then through an unheated flow developing section 
(upstream portion of the test section) which had a length greater than 40 times the test 
section’s diameter. This flow developing section had the same diameter as the heated 
portion of the test section. As the fluid flowed through the test section, the inlet and exit 
temperatures and pressures were monitored. The downstream portion of the test section 
was heated (with a heater tape) on its top half. A downstream valve was used to control 
the test section’s exit pressure. The fluid then passed through the heat exchanger where 
the energy generated in the test section was removed. For these test, the working fluid 
bypassed reservoir 2 and the charging pump and then flowed back to reservoir 1. 

Test Section 

The test section was 223.0 cm long and is shown in Figures 2 and 3. The upstream 
unheated portion of the test had smooth walls and the downstream heated portion (121.9 
cm long) had either smooth or an enhanced wall configuration (spiral fins and/or a twisted 
tape). The entire test section was insulated, and had three main ports (with a forth ex- 
tending form the center one), mounted facing downward on either end of the test section. 
These ports were used to monitor the inlet and exit fluid pressures, temperatures, and test 
section differential pressure. Each test section had twenty-eight (28) Type-K thermocou- 
ples mounted on the outside surface of the copper flow channel. Seven thermocouples were 
used to make temperature measurements ( T m ) at specific axial locations (see Figure 3) 
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on the wall of the coolant channel. Figure 3 also shows the four circumferential locations 
(^ = 0, and 7 r radians) at which wall temperature measurements were made for each 
of the seven axial locations. All measurements were made for flow conditions of: (1) 0.19 
MPa (absolute) exit pressure, (2) 22.2° C inlet temperature, and (3) a mass velocity of 
210 kg/m 2 s. These wall temperature measurements were used, along with other measured 
conditions and the data reduction analysis, to determine the unknown steady-state heat 
transfer coefficients (h, and h m ). 

Data Reduction 

The data reduction approach was based on a heated hydraulic diameter (Boyd and 
Turknett, 1989) assumption. Figure 4 shows the model used for this approach. This 
model was used to compute a circumferentially averaged heat transfer coefficient from the 
circ umf erentially averaged wall temperature. This latter temperature was computed from 
the four wall temperature measurements made on the outside of the test section at each 
of the seven axial locations. 

Briefly, this approach involves estimating the inside flow channel’s wall temperature 
by using an equivalent uniformly heated tube, whose diameter is equal to the ratio of 
four times the actual flow channel cross sectional area to the heated perimeter. This was 
done using the model in Figure 4, by accounting for the temperature drop across the flow 
channel w all , and the heat losses (convection and radiation) from the test section to the 
ambient. An iteration scheme was necessary to compute the inside wall temperature. After 
accounting for finite heat losses, the circumferentially averaged heat transfer coefficient was 
given by 
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In the above expression for the circumferentially-averaged heat transfer coefficient 
(h m ), the magnitude of the bulk temperature, T/, is dependent on T av , which is the 
circumferentially-averaged wall temperature at r = rc (see Figures 3 and 4). Tf was 
determined based on the magnitude of T av relative to the wall temperature, T WONB which 
is the temperature required for the onset of nucleate boiling. The bulk fluid temperature 
is given by, 

T _ f T f (z), for T av < Twonb'i ^2) 

^ \ Tsaty ^av — TWonb 
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For a given axial location, the measured circumferential values of T m were related to T av 
by 

T m {9 = 0) + 3 T m (9 = f ) + 3 T m (9 =*-?) + T m {9 = tt) 

Tav — g » W 


Finally, the temperature T Wonb was computed using the correlation by Frost and Dzakowic 
(Collier, 1981), 

1.0 = 8.0 StoNB P& We* 1 Ec Ja*. (4) 

Using equation (4), the superheat required for the onset of nucleation was 5.9° C (for a 
pressure of 0.19 MPa). 

In some cases, the heat transfer coefficient was not only averaged circumferentially 
but also axially. Although there were seven axial locations at which wall temperature 
measurements were made, those measurements near either end of the test section heated 
length were influenced by end losses. Hence, the averaged heat transfer coefficient, h, was 
obtained using the five central axial locations (Z2, Z3, Z4, Z 5, and Zg), 


h ^mi ~t~ 2/t m3 -f* ~h h me 

8 

so that the values of h mi ( where i = 2, 3, 4, 5, and 6) correspond to the various locations 
at Zf . 

An uncertainty analysis was developed using the above formulae to estimate the un- 
certainty, 8h m . Using the approach suggested by Moffat (1988, 1990), 8h m was found to 
be ±14.6 W/m 2 K. 
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RESULTS 


For the various internal channel configurations noted above, comparisons were made 
of the two-dimensional wall temperature distributions, axial distribution of the mean 
(circumferentially-averaged) heat transfer coefficient ( h m ), and the totally averaged (cir- 
cumferentially and axially) or overall heat transfer coefficient (h). 

OveraU Heat Transfer 

Figure 5 shows a comparison of the overall heat transfer coefficient for the four in- 
ternal configurations. These comparisons show that the spiral fins with the large pitch 
resulted in a higher heat transfer coefficient at all power levels. The discontinuities in each 
curve are due to either nucleate boiling or severe flow structure change at certain axial 
or circumferential locations for a given power level. Since the test section was horizontal 
and since the mass velocity level was relatively low (low FVoude number), stratification 
effects were expected and found to be significant. Stratification conditions reduced the 
enhancement effectiveness for all internal configurations. Preliminary estimates indicated 
that these reductions could be as high as an order of magnitude relative to vertical flow. 
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Two—D Wall TemperatMre Distribution 


Figures 6 (6a-6c) and 7 (7a-7c) show the power generation as a function of measured 
outside wall temperature at different circumferential and axial locations for the cases of 
large pitch spiral fins without and with a twisted tape, respectively. The three figures in 
each of these sets axe for three of the four circumferential locations (phi=<^ = 0, and 7r; 
see Figure 3). As <f> varied from 0 to j, the peak wall temperatures for the case of large 
pitch fins with the twisted tape were consistently higher than those for the channel with 
only the large pitch spiral fins. This was displayed more dramatically in Figure 8, which 
shows the axial distribution of Tw (= T m ) for the four circumferential locations and 
constant power. While the wall temperature distributions were essentially identical for 
^ and 7 r, there were significant differences at <f> = 0 and f (compare Figures 8a 

and 8b). It is apparent from these figures that the addition of a tape increased mixing 
and reduced stratification at downstream locations. For the case without the twisted tape 
and for <f> = 0, the wall temperature in the upstream portion of the test section between 
Zi (Z=20.32 cm) and Z 4 (Z=60.98 cm) was consistently lower than that with the twisted 
tape. However, this trend reversed downstream of Z 5 (Z=81.28 cm). Since the power level 
for the large pitch finned tube in Figure 8a (P p — 653 W) was greater than that for the 
large pitch finned tube with a twisted tape (shown in Figure 8b with Pp— 611.0 W), the 
profiles for the former case will be consistently lower than that for the latter case at the 
same power level. With the exception of the differences noted at <f> = 0 and the two 
profiles were similar. Further, for the case with the twisted tape and the large pitch fins 
and at <f> — j, the wall temperature varied axially in a periodic manner between maxima 
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of 155°C (upstream) and 90°C (downstream). Since there were only seven axial locations 
in which measurements were made, it was not possible to determine the period of this 
variation. However, the present measurements indicate that the period was less than 40.0 
cm. As implied above, the amplitude of the fluctuations decreased as Z increased. These 
latter trends were caused by: (1) periodic liquid wetting near ^ = f (off center from the 
the top of the channel) due to the swirl flow, (2) liquid entrainment into the vapor flow, 
and (3) circumferential conduction in the tube. 

FVom Figures 8a and 8b, stratification effects can be seen to have been significant 
in that: (1) the wall temperatures at <f> = 0 remained significantly above the saturation 
temperature (T sat = 43°C), and (2) those wall temperatures at <f> = 2s and tt were 
consistently below T aa t- However, the tubes with the large pitch fins (with and without 
the twisted tape) reduced stratification much more than the small pitch finned wall tube. 
As Z was increased, the stratification decreased more for the former two configurations 
than for either the smooth tube or the small pitch fins (compare Figures 8a and 8b with 
8c and 8d). Although the power levels for the latter two figures are lower than that of 
the former two, a comparison of the axial distribution of the circumferentially-averaged 
wall temperatures clearly shows significant reductions in Tw with respect to Z for the large 
pitch finned wall. In fact, the large pitch spiral fin with the twisted tape was more effective 
than all other cases. 

In all tube configurations, the wall temperature increased (from 22.2 °C) with Z near 
the entrance (Zj) and later decreased (to near 25.0°C) as Z approached Z7, near the exit 
of the test section. 
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Ctrcum/erenriafZ;/ Mean , Axially^ Distributed Heat 


Figures 8a and 8b emphasize the significance of the circumferential temperature vari- 
ations in systems with single-side heating. It is apparant that in cases where stratification 
is important, the addition of a twisted tape will, at some locations, exacerbate (compare 
Figures 8a and 8b, for <f> = j) the already large wall temperature and small heat transfer 
coefficient. However, the data also indicated that the twisted tape will enhance, rather 
that inhibit, the heat transfer in some cases. Before this enhancing effect is discussed, the 
adverse influence of the twisted tape on the heat transfer coefficient will be discussed. 

The detrimental influence on h m due to the addition of the twisted tape was empha- 
sized by considering the circumferential averaged heat transfer coefficients as a function of 
the power generation, with Z as a parameter. At each axial location, a sudden rise in h m 
was a manifestation of the inside wall temperature (computed from the measured outside 
wall temperature) exceeding the absolute wall superheat required for the for the onset of 
nucleate boiling. Relatively speaking, larger values of h m (1,500 to 2,000 W/m 2 K) were 
obtained at both the entrance ( Z \ ) and the exit (Z 7 ) than at intermediate locations. This 
is due to: (1) entrance effects, (2) the presence of the single-phase liquid at the bottom 
of the tube, (3) axial conduction losses, and (4) the absence of heating from the heater 
tape near the exit. Therefore, the data near the test section exit and entrance may not be 
representative of the actual behavior. Nevertheless, when the values of h m at the interme- 
diate axial locations were compared, one finds that the levels of h m before and after ONB 
were higher for the tube without the twisted tape. 

The above trends could possibly be reversed by reconfiguring the twisted tape. The 
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twisted tapes’s twist ratio ( tx , ratio of the axial period to the inside diameter) appears 
to be the underlying factor which could improve the enhancement capabilities of the tube 
with both fins and a twisted tape. Supporting evidence for this possibility can be observed 
by comparing either Figures 6b with 7b or 8a with 8b for <f> = As noted earlier, the 
effect of the twisted tape is to raise the local wall temperature at some axial locations and 
lower it at other locations. In cases where the wall temperature was lowered, these lower 
values (as well as the peak values) decreased with increasing Z. It would appear that if the 
period of these temperature fluctuations could be decreased, the lower levels of the wall 
temperature would prevail over a large portion of the flow channel. Lower temperatures, 
and hence larger would result due to increased mixing between the stratified fluid 
layers. This overall trend, of enhanced heat transfer accompanying reduced tj-, has been 
pointed out in the literature (e.g., see Kirishenko,1980; and Hong and Bergles, 1976) but 
has never [to the authors’ knowledge] been documented by local measurements on top- 
heated tubes. However, to verify that this is also true for stratified flows, the present work 
should be extended to include lower values of tr- 

The axial distribution of the circumferentially mean heat transfer coefficient ( h m ) is 
shown in Figure 9 for three of the four internal tube configurations. The axial trends for 
the smooth wall case in Figure 9c form a basis for the other cases. The trends in the axial 
variation of h m (Z ) become increasingly irregular as the internal enhancement progresses 
from smooth wall tube to the small pitch fins and eventually to the large pitch fins with 
the twisted tape. There was a reduction and subsequent increase in h m (Z ) as Z increased, 
which is similar to observations made by Reid et al. (1987) for the case of uniform heating. 
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FYom Figure 9c, there was a local (axial) peak in h m (Z), which moved downstream as the 
power increased. This may be representative of a slug-type flow and may be a unique 
consequence of the top-heated boundary condition. The width of the axial distribution 
and the magnitudes of h m (Z ) increased with the power. For the smooth wall case, all 
curves appeared to approach an asymptote as Z increased. The addition of enhancement 
devices disrupted these rather regular trends. The fluctuaions of h m (Z ) with increasing Z 
still existed but became more irregular and greater in amplitude. 

CONCLUSIONS 

The present work provides experimental data characterizing the localized thermal 
transport in top-heated horizontal coolant channels with enhancement devices. The present 
local wall temperature measurements form a basis for future comparisons with both three- 
dimensional numerical predictions. Such comparisons will be useful in explaining the 
underlying local flow conditions which are favorable to both local and overall heat transfer 
enhancement in top-heated configurations. This can be demonstrated in a limited way 
by noticing the behavior for h m for <f> = j in Figures 8a, and 8b. The effect of adding a 
twisted tape was to move the peak wall temperature upstream. The present cases should be 
expanded to include additional circumferential and axial resolution of the wall temperature 
variations, and comparisons at additional levels of mass velocity. The two-dimensional wall 
temperature measurements were used to determine the circumferent i ally-mean but axially- 
dependent heat transfer coefficient. The flow in the coolant channel was hydrodynamically 
developed but thermally developing with regions of: (1) single phase convection, (2) local 
subcooled boiling, and (3) a predominating stratification flow over most of the channel’s 
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length. 

The results show that the coolant channel with the large pitch spiral fins had a larger 
overall heat transfer coefficient than smooth tubes, or tubes with either small pitch spiral 
fine or a combination of large pitch spiral fins and a twisted tape. However, local measure- 
ments indicated that the effectiveness of the latter case will improve for stratified flow as 
the period of the twisted tape is reduced. Although similar observations have been made in 
the literature for non-stratified flows, the present local measurements not only documents 
t.V»ig effect, but: (1) provides a basis for comparisons with three-dimensional, two-phase, 
numerical models, and (2) forms a basis for assessing present and evolving heat transfer 
correlations. 
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Table Is Test Section Internal Configurations*. 


Tube Type 

O.D. 

I.D. 

No. Fins 

Fin Height 

Fin Width 

Fins/cm 

Spiral Fin L.P. 

1.27 cm 

0.95 cm 

16 

0.056 cm 

0.30 cm 

4 

Spiral Fin L.P. /Tape 

1.27 cm 

0.95 cm 

16 

0.056 cm 

0.3 cm 

4 

Spiral Fin S.P. 

1.27 cm 

1.13 cm 

26 

0.056 cm 

0.30 cm 

6 

Smooth Walls 

1.27 cm 

1.07 cm 

— 

— 

— 

— 


*Note: L.P., and S.P. denote large and small pitch fins, respectively. 



Figure Captions 


Figure 1. Freon- 11 Flow Loop for Both Subcooled and Saturated Flow Boiling Experiments. 
Figure 2. Cross Section of the Heated Portion of the Test Section. 

Figure 3. Wall Temperature Measurement Locations. 

Figure 4. Control Volume for the Heated Hydraulic Diameter Model. 

Figure 5. Comparison of the Overall Heat Transfer Coefficients for Circular Coolant Channels 
With Different Internal Configurations. 

Figure 6. Measured Outside Wall Temperature Axial Distribution as a Function of the Net 
Power Generation for the Spiral Fin, Large Pitch Internal Geometry at: (a) <f> = 0, 

(b) <f> = f , and (c) 4> = 7r. 

Figure 7. Measured Outside Wall Temperature Axial Distribution as a Function of the Net 
Power Generation for the Spiral Fin, Large Pitch With a Twisted Tape Internal Ge- 
ometry at: (a) ^ = 0, (b) <f> = j , and (c) <f> = 7r. 

Figure 8. Axial Distribution of the Wall Temperature for the Four Circumferential Locations 
for the: (a) Large Pitch Spiral Fins, (b) Large Pitch Spiral Fins With Twisted Tape, 

(c) Small Pitch Spiral Fins, and (d) Smooth Tube. 

Figure 9 Axial Distribution of the Circumferentially Averaged Heat Transfer Coefficient for. 

(a) Large Pitch Spiral Fins With a Twisted Tape, (b) Small Pitch Spiral Fins, and 
(c) Smooth Wall. 
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Figure Is Freon-11 Flow Loop for Both Subcooled 
and Saturated Flow Boiling Experiments. 





Thermocouples 
(4 locations) 
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Figure 2- x press Section of the Heated Portion of the Test 
Section. 
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Figure 3.t Locations in which Local Outside Wall Temperature 
Measurements were Made. 
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Figure 8. (b) Large Pitch Spiral Fine With Twisted Tape, 
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Figure 8. (d) Smooth Tube. 



TOP-HEATED 1.27cm SPIRAL FIN L.P. W/TAPE; HZ 1122 
a = 210 ka/m's: INLET TEMP. = 22.2 C 
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Figure 9 (a) Urge Pitch Spiral F.ns With a Twirled Tape, 
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APPENDIX: Uncertainty Analysis 


This study has resulted in the determination of circumferentially mean and axially 
(and circumferentially) mean heat transfer coefficients from measurements of: (1) the test 
section outside, local wall temperature (2) the freon-11 flow conditions which include flow 
rates, exist pressure, and inlet and exit bulk temperatures, and (3) the ambient tempera- 
ture. The relationships between these quantities are summarized in equations (1) through 
(5). Following the approach outlined by Moffat [7,8], the uncertainty in the heat transfer 
coefficient, h m , is 

Sh m = [(|5=. SXi) 2 + (~6C,) 2 i ,/2 (A. 1) 

where Xi represents all independent variables, and C< represents corrections used to ac- 
count for calibration defects, system-senson interactions and system disturbance errors. 
The double indicies in eq.(A.l) imply summation over all independent or correction vari- 
ables. If effects of C{ are neglected, 6h m is given by, 

«„*[( iAf + { ^«B) 2 ] ,/2 (A2) 

where A and B are given in equation (1). Since A and B are not actually independent 
variables, thier relationship with these variables are as follows: 

6 A = Function(T m , PpiToo), so that (-4.3) 


f)A 2 8A 


) 2 + (^^ £T °° ) 2 ] 1/2;an< * 
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SB = Function(T 00 ,T m ,Tf,P p ),so that 


U-4) 


SB = iMtP,) + (^ST m ) + (%-&<*? + ] 1/2 


dB 


dB 


dB 


dP v 


dT, 


dT a 


dT f 


FVom equation (1) the partial derivatives shown in equation (A. 3) are given by: 


dA 

dP p 



r_c 

A s t a 


)[ln& 

r c 


+ 


kc 

hoo r D 


], 


(A.5) 


dA _ i dhpo _ kc 

dfZ~~r A hJr D dT m r A 


dA _ _ k^kcBi dhpp k c 
dToo dToo t a 


In the above expression, As = 2 tt DL, where D = 4 r A . The partial derivatives involv- 
ing B (see equation A.4) can be obtained using equation (1); i.e. 


dB 

dP p 


r_c 

A 3 kB 


[In(jr-) + 

rB k A r A 


(A.6) 


dB 

dT m 


= $ - (T m - Too)( 


kc x dhoo 
hlo ro dTm 


+ ( T f - T^) 


kc dh^o 
k 2 oo r D dT m 


; where. 
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* = (In(^) + 

r c 


kc 

h<x> I'D 


+ 


k c 

kB 



dB 

dT^ 


-$ + /n(— ) + 

r c h'oci'D 


+ (T, - Too) 


k c dh c 
hoo^D dT c 


and 


dB 

dT f 



kg 

hoo r D 


] 


The local fluid bulk temperature was estimated via the First Law of Thermodynamics, 
the measured values of mass velocity, net power generation, inlet temperature, pressure, 
and the axial location;i.e., 


Tj = Function(G,Pp,Ti n i e t,P, Z) 


(A. 7) 


so that, 


stM(§-6p p ? + 


where 


dT_ '/ _ ( 0, for T w < T WO nb] 
dT sat I 1> f° r Pw ^ Pwonb'i snd 
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ZPr, 


Tf - TinUt + TLr t GC 


A P 


The above partial derevatives axe given as, 


dTf _ Z 
dP p 7r Lr\GCp 


dT inUt ’ 


(A.8) 


dT f P P , 

dZ LA c GCp ’ 


and 


5T/ __ ZP p 
dG ~ 7T Lr\G 2 Cp 

A complete energy balance of the test section was performed. This required appro- 
priate definition for /loo- The effective external heat transfer coefficient, hoo, was defined 
in terms of the external natural convective heat transfer coefficient (h c ), and the external 
radiative heat transer coefficient (h r ) by 


hoo = h c + h r ^h c + <r(T 2 + Too 2 )(T a + Tool ( A . 9 ) 

where a is the Stafan Boltzmann constant and T t is the outside surface temperature of 
the test section. The temperature, T a , is related to T m by the approximate relationship, 
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(.4.10) 


2itL(T m - T.) 2nL(T m - T„) 

**•(£) 

Equation (A.10) can be used along with (A.9) to solve for T„ iteratively for each 
measured value of T m . An itenative analysis is necessary because is a nonlinear function 
of T a via h r . The natural convective heat transfer coefficient, h c , was determined from the 
following correlation (for horizontal cylinders) 


h e 


2r D 


where 


(All) 


Nud — {0-6 -f 0.38i?a£)^[l + ( p^ m ) TC ] ^} 2 > 


where 


Rap = 


PgjT, - Too )!) 3 

I/ 2 


The partial derivatives of h 0 0 with respect to T and T m appear in equation (A.6) and 
are given by 


dh op 

dT^ 


p- + c[2 (r.|A + x.xr. + r.) + (T . 2 + t to 2 x ^ + i)i 


(A12) 
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dh op 

dT m 


dh c 

dT m 


+ cr[3T, 2 + 2 T'Too + Too 2 ] 


dT, 
dT m ’ 


where and were obtained from equation (A. 10). The resulting expression for 


ar c 


lfc is 


ar. 


*cV 

*»<g) 


b 5 - 




fp*c 


¥ 


dT m 


ar "> + ''P 1 '- 2 + 2T - r ~ + r ~ 2 i 


U-13) 


where 



1 

^oo I'D 


Sample calculations indicated that is about 1.9 x 10 2 at higher power levels 
(P p = 657.5W), and should be much smaller (nearly zero) at lower power levels. An 
expression similar to equation (A. 13) can be obtained for ~§y^ from equation (A. 10) and 
(A. 12). Sample calculations have shown that h 0 0 is nearly constant for small changes in 
Too- Assuming this to be the case, especially at higher power levels, 


dT, 

dToo 


In^P- 

IQ 

k r B$ 


(A.U) 


which is slightly less than unity. 

One of the last variables for which an uncertainty must be estimated is the net power 
generation. Since T„ was near Too , <5P p will be strongly dependent on the measured 
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electrical current (I) and voltage (V); i.e., 

IP, = [(|f«) 2 + (|^0 2 ] 1/2 (^-15) 

Here we assume that the uncertainty for the partial derivatives of the power and total 
power generation are identical. Noting that Pp) tota i = IV, 

6P p = [ V 2 61 2 + I 2 SV*] 1/2 (A.16) 


Sample Calculation 


Since the primary objective of this work was to obtain heat transfer coefficients for 
enhanced and smooth surfaces in the boiling regime, one power level which resulted in 
localized boiling has been considered. The above equations are being used to develop a 
computer program so that the uncertainty at every power level can be evaluated. An 
undergraduate mechanical engineering senior was employed for this purpose. 

In particular, selected conditions at a typical high power level was considered for sam- 
ple calculations. The test section with large pitch fins for <f> = 0, Z = Z4, and P p = 657.45 IT 
[see reference A-l, p. 232]. In addition to G = 210 kg/m 2 s, Z4 = 60.96 cm, Tf = 314.8 A , 
Too = 295 K, and T m = 376.87 K. 

In the original data reduction used in Part I and [A-l], T t was assumed to be 7 ’ m + 7 °° 
and was used in equation (A. 9) to compute h T . The heat transfer coefficient, h c , was 
assumed to be 4.24 W/m 2 K for that data reduction. However, the present computa- 
tions show that T, should be 296.5 K [obtained from equation (A. 10)], h c should be 
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28.6 Wfm 2 K, and should be approximately 32.3 W/m 2 K. This change has a sec- 


ondary effect on the uncertainty analysis since was small compared to which 

is a controlling parameter in this analysis. 

The following is representative of the specific characteristics of the test section dimen- 
sions, fluids, and thermophysical properties: r A = 6.033 x 10 -3 m; t b = 6.35 x 10 3 m; 
r c = 6.98 x 10 -3 m; r D = 3.2385 x 10“ 2 m; D = 4r^; L = 1.219 m; k A = 382.7 W/mK] 
k B = 1.038 W/mK; k c = 0.0303 W/mK ; k ^ = 31.3 x 10~ 3 W/mK (Air); i/ TO = 
0.2945 xlO -6 m 2 /s; (Air); /? = 3.3 xlO -3 K~ l (Air); Pr ^ = 0.704 (Air); p « 1,479 kg/m 3 
(freon-11 density); C p « 0.89 kJ/kgK (freon-11 ). 

To complete the specifications for the sample calculations, the measured uncertainties 
of the independent variables must be determined. The uncertainties for Z and G were 
: SZ = 0.0016 m, and SG = 8.4 kg/m 2 s , respectively. The uncertainties in measuring 
all temperatures were assumed equal and will be denoted by ST. If SP p = 0.1 W , ST = 
10~ 2 K, then 6 A = 762.3 W/m 2 , SB = 0.027 K, = 9.92 x 10' 3 K~\ = 

—462.6 W/m 2 K 2 . This results in 


6h m w ±14.6 W/m 2 K. 


(A.17) 
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Discussion 


The value of Sh m shown in equation (A. 17) is strongly dependent on ST and 8P p . This 
value would be much smaller for lower values of ST and SP p , or significantly larger for higher 
values of these quantities. In the expression for SB , the quantities which have secondary 
influence were: |^- = 2.8 x 10 -2 K/W] and = -1.02 x 10 -3 . Those quantities 
with significant influence were: = 1.534, B = 100.8 K, and = —1.549. There 

is great concern in the large magnitude of this last quantity; i.e., • It’s magnitude 

arises from many influneces associated with the magnitude of Tj (the local bulk fluid 
temperature). For T/ < T aai , ^p- — 0 and the uncertainty in T/ is ST as noted above. 
However when T f > T tai , there is a possibility that ST f may become large. For example, 
15^ « 180.9 °C/MPa « ^4- for boiling around the entire channel. If this condition 
occurred (e.g. , with uniform heating) with only a moderate uncertainty in P (i.e., SP — 
0.014 MPa), STf could be as high as 2.4 K. However, in most experiments the bulk fluid 
temperature remained below the saturation temperature even when there was localized 
boiling. Under these circumstances, STf would be slightly greater than ST, provided 
SP < 6.8 x 10" 5 MPa. 

The un certainty, Sh m , is also dependent on the parmeter A, which has an unusually 
large magnitude of 4.7 x 10 6 W/m 2 for sample case considered. The uncertainty, SA, is 
principally affected by P p ; so, the effects of -gfc = = 116.0 W/m 2 K are secondary. 

Accordingly, = 7,623 m -2 . 
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Experimental Design 


Many improvements for the existing experimental loop have been enumerated in [A-l]. 
In addition to those, the heat exchanger for controlling the test section inlet temperature 
was relocated upstream of the test section. Additional verification runs are planned in 
anticipation of additional funding to study the vertical flow geometry which has been 
proposed recently [A-2]. 

In addition to the above discussed uncertainty sensitivity, this uncertainty analysis 
suggests some subtle changes in the test section which will reduce the uncertainty in 8h m . 
In addition to requiring low £T m and k c , if r D /r c was as large as possible, both and 
would be significantly reduced. This occurs because both of the latter quantities 
are inversely proportional to B and B ^ , respectively. The quantity, B, would increase 
further if T tat could be selected to be closer to T^. The latter requirement would have a 
significant effect on reducing A, and hence and 8h m . Another very large contribution 
to 8h m was due to 8 A. The quantity most significantly affecting 8 A was ^pn. This latter 
quantity can be significantly reduced by choosing a low thermal conductivity material 
for the test section (i.e., k A ). For example, if k A is reduced by an order of magnitude, 
8 A will be reduced similarly (i.e., by and order of magniftude). Any one of these effects 
separately will reduce 8h m by a relatively small amount. However, by employing all of 
them, significant reductions will occur. 
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ABSTRACT 

A heat transfer model for a single-side, heated, circular coolant channel has been developed 
and can be used to determine accurate mean heat transfer coefficients from measured wall 
temperatures. In addition, the circumferential heat flux distribution on the inside wall of 
the coolant channel was obtained. In previous experiments, seady-state wall temperature 
measurements were made at various circumferential locations on the outside surface of the 
channel which was heated from one side and insulated on the other side. The channel was 
assumed to be cooled internally by a fluid flow field with a bulk temperature, TJ,. In the 
present work, cases of uniform and nonuniform heat flux were considered for the heated 
side. The results were persented for the Biot number varying from 10 -3 to 10 3 , and aspect 
ratio ( R 0 , outside to inside radii) varying from 1.04 to 3.0. Using the formulation, the 
circumferentially mean heat transfer coefficient (h m ) can be obtained from specified or 
measured values of the outside surface temperature T(R 0 , 4 > ')i &nd T&. If measurements of 
the local fluid film temperature are made available, the present formulation can be used 

1 Honeywell Endowed Professor of Engineering 

2 Graduate Student 
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to obtain the circumferential variations of the heat transfer coefficient, 

Since the results are closed-form analytical expressions, they can be used conveniently 
for data reduction or for boundary conditions for numerical analyses of the flow field. In 
cases where axial (local) wall temperature measurements are made along the length of the 
channel, the present results can be used locally to approximate variations of the mean heat 
transfer coefficients in the axial direction. 


NOMENCLATURE 

a„ Integration parameter; see equation (4) 

Bi Biot number, hm r l / k 

Bi n (Bi — n)/(Bi + n) 

b n Integration parameter; see equation (4) 

h Local heat transfer coefficient 

h m Mean circumferential heat transfer coefficient 

k Wall thermal conductivity 

m n =(-h m rpnfcrf n-1 )/(-h m r 1 n + fcnrj* -1 ) 

n Eigenvalue 

q ” Circumferentially varying, outside surface heat flux 

q 0 Constant heat flux 

r Radial coordinate 

R r/r a 

R a Aspect ratio, r 0 /ri 
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r 

r 0 


r i 
T 
Tj 
T* 


Radial coordinate 
Outside radius 
Inside radius 
Local wall temperature 
Bulk fluid temperature 

Dimensionless temperature (referred to as T STAR in plots), B/(q 0 r 0 /k ) 

Average circumferential value of T * (referred to as AVERAGE T STAR in plots) 


Greek 

B(r,<f>) Temperature excess above the bulk fluid temperature, [T(r,</>)-T6] 
<f> Circumferential coordinate directed counter-clockwise 1 


1 The origin for <f> is at: (1) the top of the vertical plane (see Figures 2 and 3) for all 
plots, but at (2) the right-side of the horizontal plane for equations (1) through (14). 
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INTRODUCTION 


Recently, much attention has been given recently to single-side heated coolant chan- 
nels. Most practical applications require forced convective cooling (single-phase or two 
phase) and have coolant channels which have drastic circumferential variation in the ap- 
plied heat flux. Sometimes, the axial heat flux variations are secondary compared to the 
circumferential variations. In cases where this is not true, the present data reduction tech- 
nique must be applied at each local axial location. Because of the simplicity involved, 
many previous investigators and designers have approximated single side heated coolant 
channels with an equivalent uniformly heated channel. Since this approximation may not 
apply in all cases, more investigators are now interested in including the effect of single-side 
heating in their investigations. W ith the addition of this needed complexity to man\ ex- 
periments, new and convenient approaches must be identified to reduce heat transfer data. 
This need is not restricted to experimental investigators. Many numerical analyst neglect 
the effects of the congugate nature of their analysis by approximating heat flux boundary 
conditions rather than using the actual boundary conditions which may have been com- 
plicated by multidimensional conduction (and/or radiation) in the solid substrate. Within 
the scope of the present analytical formulation, the analysis which follows provides a means 
for those analyst to include more realistic boundary conditions for their study of the fluid 
flow through non-uniformly heated circular channels. 

In many engineering designs, various heat transfer enhancement devices (e.g. fins, 
twisted tapes, combinations of fins and twisted tapes, etc.) cause very complicated internal 
velocity and temperature distributions. Consider a circular coolant channel which is: (1) 
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heated from the outside via a nonuniform heat flux, and (2) cooled from the inside by a 
complicated convective flow distribution at a bulk fluid temperature, T h . Such convective 
flows are often times too complicated to analyze numerically, especially if the flow is 
turbulent and/or involves a phase change. The above described configuration models, in 
two-dimensions, nonuniform heat flux configuration found in electronic devices, reactor 
components, space cold plates, and many other energy management systems. 

For most thermal engineering designs, knowledge of the circumferentially mean heat 
transfer coefficient, h m , is essential and knowledge of the local circumferential heat transfer 
coefficient, h (<£), and hence the local wall temperature would be desirable. Although 
previous models [1-4] have been developed for the case of nonuniform surface heat flux, 
they have either not adequately accounted for the nonuniformity of the heat flux or they 
were based on laminar flow. Apriori assumptions of such distributions are usually not 
appropriate. 

The objectives of this work were to: (1) establish a convenient technique for deter- 

mining h m from measured or set values of the outside channel wall temperature, T(R 0 , d), 
and (2) propose a technique for using the external heat flux distribution to infer the inside 
wall heat flux distribution. These set values for T(R 0 ,<f > ) or q ” may either be experimental 
measurements or critical design constraints. 
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PROBLEM SPECIFICATION AND FORMULATION 


The model to be studied consists of a nonuniformly heated coolant channel (see Figure 
2) which is heated from one side and insulated (perfectly) on the opposite side. The ther- 
mophysical properties of the channel material were assumed constant. The mean circum- 
ferential heat transfer coefficient, h m , was assumed to adequately describe the convective 
heat transfer from the inside surface of the channel to the bulk fluid. Although not used 
directly in this formulation, this does allows for the existance of circumferentially varying 
heat transfer coefficient, This latter coefficient can only be found by determining the 

circumferential variation in the ’’local” fluid film temperature necessary for the definition 


of h(4>). The wall temperature is charaterized by the following differential equation: 


1 d 2 T ld_ ( dT, = 

r d4 > 2 r dr dr ’ 


However, note that for the formulation, the circumferential coordinate, <j>, is directed 
counter-clockwise with its origin ( <j> = 0) on the right side of a horizontal plane. The 
origin for <f> will be changed to that shown in Figures 2 and 3 for all plots in Figures 10 
through 13. The boundary conditions axe 


T(r,4>) = r(r,<£ + 27r), 


ar(r, 4 ) _ ar(r,<i + 27r) 


-h rn [T(r,<t>)-T b ) + kj-(r,<i>)=0, and (2b) 

dT(r 0 ,<f>) _ { —q n ((f>), 0 < <f> < tt] (0 , 

- k dr" -\ 0 , ( ) 

After these equations were nondimensionalized and restructured using the standard sepa- 
ration of variables technique, the following temperature profile was obtained: 

OO 

0 = 6 0 + ^(r“ n - m n r n )(a n cosn<f> + b n sinn<f>), and (3) 
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6o = A 0 [ + InR ], 

where n=l,2,3,4,5 , which reflects all boundary conditions having been applied except 

the last. The first term, 6 0 , would have described the radial temperature profile for the case 
of a completely uniform heat flux around the circumference of the flow channel. The second 
term is therefore the required correction due to the nonuniform heat flux. Application of 
the last boundary condition results in expressions for the constants a n and b n in terms of 
the nonuniform heat flux distribution; i.e., 


(r 0 n 1 + m n r" 1 )a n = — f p-cosn(f>d4 > , 

7rn J 0 k 


(rj n_1 + m n r n 0 1 )b n = — [ and 

Jq K 


( 4 ) 


r 2n 


7T q " (m 
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Case I : q” (<f>) = q 0 


The case of q (4>) having a constant magnitude of q 0 on one side of the channel and a 
zero heat flux on the remaining side is applicable to electronic cooling, space cold plates, 
and is an approximation for reactor component heat loading and many other applications. 
For the latter application, a sinusoidal heat flux distribution, which will be covered in Case 
II (below), is more applicable ot fusion reactor components. For the case of a uniform heat 
flux on one side of a channel, the dimensionless temperature distribution is: 


rp*( ft M 6 JSjS. _ If 1 I 1 „ m | V 2sinn(fr -R” [1 - Bi n R 2n ] ( ) 

T 2 \ Bi + ^ + v n 2 Rn [i + Bi n R7 2n ]' 

where a n =0, A 0 - and b n = n=l,3,5,7.... It is instructive to examine the 

validity of this equation with respect to the limits of the Biot number, Bi. For R=l, 

Asinn<f> 




( 6 ) 


2Bi 7 mR%[(Bi + n) + (Bi — n)R 0 2n ] 

As Bi approaches infinity, the convective transport mode predominates. In that case, both 
terms in equation (6) approach zero and this implies that T(r, <^) approaches Tf > . In a 
similar fashion, T(r 0 ,<f > ) approaches a temperature above Tj as B, approaches infinity. 
The radial temperature gradient is given by 


dT* _1_ ^ 2 sinnt iT -1 {{Bi + n) + (Bi - n)R~ 2n ; 

-\R,<p) — / ,, z>n W r>: i „\ i / D. „\D- 2nl 


(7) 


2R ' ^ nn R? [(Bi + n) + (Bi — n)Ro) 

The temperature gradients at R=1 are especially important since they represent the ef- 
fective circumferential heat flux at the inside wall of the coolant channel. Therefore, for 


R=l, 


dT* 1 4 sinrufi Bi 

dR 2 ttuI?” [(Bi + n) + (Bi — n)Ro 2n ] 
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For Bi=0, the lowest achieveable dimensionless temperature gradient at the boundaries 

occurs under the conditions of very high channel thermal conductivity or very low heat 

transfer coefficient. Therefore, is given by 

+ forR = R c and any B,; 


dT‘ 

dR 


l 

2 » 


for R = 1 and for Bj = 0; (8) 

I + V°° , ABisin n ± __ f or R = l and any Bi. 

The first equation in equation (8) should satisfy equation (2c). 

Two additional relationships will be of interest: (1) the average circumferential di- 
mensionless temperature excess, T* V (R), and (2) the ratio of the inside wall temperature 
to the outside wall temperature. T* v is defined as 


r 3ir 

TU(R) = - / or (9) 

7T 

T..(R) = \{j- + InR] (10) 

Therefore, if Tav is measured experimentally the mean heat transfer coefficient, h m can 
be obtained from equation (10). In cases where measurements of T(R 0 ,<f > ) are made, it 
may be important to relate this measurement directly to T(l,^). In this instance, it will 
be important to know the ratio of the inside to outside channel wall temperatures. These 
quantities are related to the ratio of T*(l,<f>)/T*(R 0 ,<i>)- It is very important to note that 
this ratio is completely independent of the magnitude of q c " and is only implicitly related 
to the thermal conductivity of the channel, h m , and the geometry via the Biot number; 


i.e., 


I . oc 4sinn4> 

T*(l, <ft) m + 2^n=l 


T*( 7? A>\ 1 r 1 . j t-> 1 i v 400 2 sinn<t> (Bi — n)R„ 2n ] 

1 {Ko,<P) 1[ + InRo] + Sn-l 

As will be displayed shortly, the above equations can be used to either: (1) obtain the 
temperature distribution for given values of q Q , h m , k, geometry, etc., or (2) obtain the 
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mean circumferential heat transfer coefficient, h m , from measurements of the outside wall 
temperature for given values of q 0 , k, and the geometry. 


Case LL 1 q ” (<f>)=q 0 sin(<f>) 


There are many applications where the flow channel is irradiated from one side by a 
heat source, and shielded on the other side. Equations (3) and (4), which apply equally 
to this case, reduce to the following after equation (2c) is applied by allowing for q” (<f>) = 
q 0 sin<f>, 


T-(RA) = l{lj-+ln R \+ljr 


n= 2 


(R-^Bin - j? n ) {(— l)"- 1 -!} 
2 n(R7 n Bi n - R%y n- 1 


+ 


{(-l)n-H _ 1} 
n + 1 


cosn<$>\ 



]sinn<f>}, 


(12 a) 


5T* „ 1,1 1 ^ (R- n - 1 Bi n + R n ~ 1 ) t 1 1 , , 

dR ( R ' ~ 7T { R + 2 ^ ( Ro n Bi n - R ?) l (n — 1) + n + ^ cosn<f> 


* r R - 2 • 

2 l /e- 1 + j? 0 {|g} 




(126) 
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Formulation Verification 


The heat tranafer coefficient h m can be defined as, 


which reduces to, 


hm — 




1.0 = 


[T(M)-T fc ]’ 

•B . -1 ^(M 1 


or 


T*(M) 


(13a) 


(136) 


1 , v^°° 4BtJinn^ 

2 + 2^n=l 7rnfig[(g i + w)+(B. -n)B7 l; *) 


1 = 


(14) 


1 . 4Bisinn<f> 

2 + ^>n=l [(g, +n )+(B.' -n)Rr z ") 

This equation represents one of the boundary conditions. All computations must exactly 
satisfy this equation for all values of , R 0 , and tf>. This is just one of many ways the 
computed results were verified. 


RESULTS* 

The results axe presented in a dimensionless form involving parametric variations of 
the averaged wall temperature, wall temperature distributions, and the radial tempera- 
ture gradient at the inside surface of the coolant channel. This gradient is most important 
because it is directly related to the local (circumferential) heat flux. For cases involv- 
ing a non- uni form external heat flux distribution, most designers are quite interested in 
the corresponding non-uniform heat flux distribution at the inside boundary. This heat 
flux is equal to —k 9T Q'*K A closed-form solution for the temperature gradient in this 
* See the nomenclature list to note the differences in the origin for <j> for the 
formulation equations and the plots, which are contained in this section. 
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latter equation has been given above, and corresponding circumferential distributions are 
presented below as functions of Bi and R 0 for the case of q = q 0 - 

Pa rametric Temperatu re Distribution 

Closed-form two-dimensional analytical temperature distributions have been obtained 
for two cases involving single-sided heat flux distribution. In both cases the heat flux was 
applied on one sidej i.e., half of the perimeter of the circular coolant channel was heated 
while the remaining half was insulated. In one case the heat flux distribution was constant 
(see equation (5)); and in the other, a sinusoidal heat flux distribution was used with the 
maximum heat flux occurring on the plane of symmetry (see equation (12) and Figure 2). 
In the latter case, a zero heat flux existed from <j> - f to <f> = This analytical approach 
is an improvement over the previous data reduction technique, which employed the heated 
hydraulic diameter model used in Part I. 

The improved technique for determining the circumferentially mean heat transfer co- 
efficient, h m , is based on equation (10). This simple equation relates the circumferentially 
averaged temperature at any R to Bi, and hence h m . In the experiments described in Part 
I, local temperature measurements were made at R = R 0 , and the averaged temperature 
was computed from these measurements. Figure 10 shows the dimensionless, averaged 
temperature as a function of both Bi and R a . For the experiments in Part I, Bi was less 
than 3.0 x 10 -3 and R 0 was near 1.1. From Figure 10, this indicates that measured values 
of T* v above 0.5 give good resolution for corresponding values of B, (or h m ). 

The local variations of the coolant channel wall temperature are presented in Figures 
11 (<f> = 0 and 7 r) and 12. Figures lla-c (4> = 0; center of the heated side) appear to show 
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that there are negligible variations of T* for all R and Bi. However, closer inspection is 
provided by Figures lld-f, where it is shown that at all values of R 0 < 1.04 and Bi < 0.01, 
there are negligible variations of T* with respect to R. Figures llg-1 show similar results for 
<f> = n (center of the insulated side). The circumferential variation of T* is shown in Figure 
12 in terms of the ratio of T*(R = 1, <f>)/T*(R = R 0 ,<f>) for R 0 = 3.0. The maximum value 
of this ratio is 1.0, which is approached when: (1) Bi becomes small (i.e. < 0.005), and 
(2) as R 0 approaches 1.0. The minimum value of this ratio always occur at <f> = 0 (Heated 
side). This corresponds with T(R 0 ,<f > ) being a maximum and T(l, 4>) being a minimum. 
The maximum value of the ratio occurs at <f> = 7r, which is on the line of symmetry of the 
insulated portion of the channel. 

Inside Wall H_eat Flux Distribution 

The heat flux distribution at the inside wall of the coolant channel is proportional 
to the local radial temperature gradient at R=1.0. Figures 13a-c show that there are 
no substantial heat flux variations on the inside wall when: (1) Bi < 10 -3 forR 0 = 

1.04, (2) Bi < 10" 2 for R 0 = 1.34, and (3) Bi < 5.0 x 10" 2 for R 0 = 3.0. Notice 
that this Bi limit increases with R 0 . Hence for Bi greater than the above limits, the 
circ umf erential variation of the inside wall heat flux becomes progressively larger as Bi 
increases. 

Because the analytical results are in form of infinite series, the convergence of the 
series requires more and more terms: (1) as R 0 approach 1.0, and (2) as B t becomes larger 
(e.g., greater than 0.2 for i? o =1.04, and greater than 10.0 for i? 0 =1.34). Let e be defined as 
the difference in the infinite series of equation (8) which is computed by using ”n” -terms in 
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one case and ”n+l” terms in the other. When e = 10 -3 , poor convergence resulted for the 
above noted intervals of £,. The results in Figures 13a-c required e = 10“ 7 for acceptable 
convergence. A number of approaches were used to check and verify the results. As noted 
above, convergence was assumed when all boundary conditions were computed within a 
minimal error, the boundary condition shown in equation (2b). 


CONCLUSIONS 

An improved technique has been developed for obtaining circumferentially mean heat 
transfer coefficients from measured wall temperature data for circular coolant channels 
which are heated extermally from one side and insulated externally on the other. Although 
two cases for the outside heat flux distribution were considered (constant and sinusoidal) 
a parametric analysis was presented for the cases of constant heat flux. For both cases, 
closed-formed analytical solutions were obtained. 

Detail results were presented' for the case of constant outside heat flux. The technique 
has good resolution for B t < 0.5 and should result in minimal error in the heat transfer 
cefficient for a dimensionless temperature, T*, greater than 0.5 for R a < 1.5. Quantitative 
profiles were computed for the effective inside heat flux circumferential distribution, which 
is proportional to ■ Significant circumferential heat flux variations occurred at 

higher Bi as R 0 increased. For example, significant variations occurred for: (1) R 0 = 1.04, 
when Bi > 10 -3 ; (2) R 0 = 1.34, when Bi > 10 -2 ; and (3) R 0 = 3.0, when Bi > 5.0 x 10~ 2 . 
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Figure Captions 


Figure 10. Circumferentially Averaged Temperature, at R = R 0 Versus the Biot 
Number, Bi. 

Figure 11. Local Dimensionless Temperature Versus R With Bi as a Parameter for: 

(a) # o =1.04, and <f>=0 (center of the heated side); (b) #0=1.34, and 
4>=0] (c) #o=3.0, and <f>= 0; (d) # o =1.04, <f>= 0, and reduced range for T*; 
(e) #o=1.34, 4>= 0, and reduced range for T*; (f) # o =3.0, <j>= 0, 
and reduced range for T*; (g) # o =1.04, <f>=ir (center of the insulated side); 
(h) #o=1.34, and 4>=ir, (i) # o =3.0, and <f>=ir, (j) # o =1.04, 
and reduced range for T*; (k) # 0 = 1.34, <£=7r, and reduced range for T*; 
and (1) # 0 = 3.0, <^=7r, and reduced range for T*. 

Figure 12. Circumferential Variation of the Ratio of the Inside to Outside Wall 
Temperatures as a Function of Bi for # o =3.0. 

Figure 13. Circumferential Variation of the Local Dimensionless Radial 

Temperature Gradient as a Function of Bi for: (a) # o =1.04, (b) # 0 =1.34, 
and (c) #o=3.0. 
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